Constraints on the primordial curvature perturbation from primordial black holes 
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We calculate the constraints on the primordial curvature perturbation at the end of inflation from 
the present day abundance of Primordial Black Holes (PBHs) , as a function of the reheat temper- 
ature Trh. We first extend recent work on the formation of PBHs on scales which remain within 
the horizon during inflation and calculate the resulting constraints on the curvature perturbation. 
We then evaluate the constraint from PBHs that form, more conventionally, from super-horizon 
perturbations. The constraints apply for Trh < 10 s GeV and the inclusion of sub- horizon PBHs 
leads to a limit which is roughly three times tighter than the bound from super-horizon PBHs. 

PACS numbers: 98.80.Cq YITP-06-66, |astro-ph/0612379| 
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The spectrum of the primordial curvature perturbation is accurately measured on large scales by observations of 
the Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) surveys. However, on small scales 1 the 
r T ^ only constraints come from avoiding the overproduction of Primordial Black Holes (PBHs). 

Typically PBHs are light black holes which may form in the early Universe, in particular via the collapse of large 
density perturbations [ITS]. If t ne density perturbation at horizon entry in a given region exceeds a threshold value 
(which is of order unity) then the region will collapse to form a PBH with mass roughly equal to the horizon mass. 
\ There are tight limits on the initial abundance of PBHs from the consequences of their expected evaporation via 
Hawking radiation and their present day abundance [!, U, [H, 0, 0, H, OH EH El, EH (for a review of these constraints 
^\ ' see Ref. |14j). These limits can in turn be used to constrain the amplitude of the primordial perturbations on small 
r*"** . scales and hence constrain models of inflation over a far wider range of scales than those constrained by CMB and 
LSS observations [H, El, E3, El, ES H3| • To date these calculations have been carried out in terms of the density 
perturbation. In this paper we calculate, for the first time, the constraints in terms of the amplitude of the primordial 
• curvature perturbation power spectrum. 

The standard PBH formation calculation applies to scales which are well outside the horizon at the end of inflation, 
with PBH formation subsequently occurring on a given scale, if the fluctuations are sufficiently large, when that 
scale re-enters the horizon after inflation. Lyth et al. [2l| extended the calculation to smaller scales on which the 
perturbations do not exit the horizon during inflation, and therefore never become classical. We use the sub-horizon 
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PBH formation calculation of Ref. [2l| to calculate the constraint on the primordial curvature perturbation from the 
present day abundance of sub-horizon PBHs and compare this with the conventional constraint from scales which exit 
^ ■ the horizon during inflation. 

Throughout this work we use 'super-horizon' ('sub-horizon') to refer to scales which do (do not) exit the horizon 
during inflation. Note, however, that in the super-horizon case PBH formation occurs when the scale re-enters the 
horizon after inflation. In Sec. [n] we outline the calculation of the evolution of the Bardeen potential after the end of 
inflation during radiation domination from the initial conditions set during the inflationary epoch, and the formation 
of PBHs on sub-horizon scales. In Sec. lIIII we calculate the constraints on the primordial curvature perturbation from 
the present day density of both super- and sub-horizon PBHs before finally discussing our results in Sec. II VI 

II. PBH FORMATION ON SUB-HORIZON SCALES 

During inflation the vacuum fluctuations of any light scalar field are promoted to a classical perturbation after 
horizon exit [H, [23[ . The quantum to classical transition of the fluctuations can be described using the Heisenberg 
picture. In the Heisenberg picture the mode function of the field perturbations, 54>k-, becomes real up to a constant 



Here and in the following by "small scales" we mean scales that are smaller than the horizon size at decoupling. 



phase transformation well after horizon crossing (see Eq. (|2.4p below), and therefore the time evolution of the cor- 
responding operator S$k(t) = 5(f>k(t)au + 5(f>t(t)a_ k , where fik and are the annihilation and creation operators 
respectively, is trivial ensuring that it has a steady eigenstate over subsequent times. For fluctuations on scales that 
never exit the horizon during inflation, the evolution of the field operators, 64>k(t), can not be regarded as classical 
because in this case the field operators do not commute at different times. However, if the perturbations are suf- 
ficiently large, then PBH formation could still occur at the end of inflation on sub-horizon scales if we regard the 
formation of a black hole itself as a measurement of the overdense region. In Ref. [2lj | it was shown that the formation 
of sub-horizon PBHs in this context should be formulated in terms of the Bardeen potential $, rather than using, 
for example, the density contrast. To find the abundance of sub-horizon PBHs we therefore need to calculate the 
Bardeen potential immediately after inflation, for scales which remain sub-horizon during inflation. This was done in 
Ref. 21f. We summarise here the essential details of the calculation. 



A. Evolution of the Bardeen potential 

We write the metric with scalar perturbations in the longitudinal (or Newtonian) gauge as [24[ 

ds 2 = a 2 (t) [- (1 + 2*) dr 2 + (1 + 2$) 5 tj dx l dx j ] . (2.1) 

In this gauge the metric perturbations coincide with the gauge invariant Bardeen potentials defined in (25j , W = $>a 
and <f> = For a scalar field the anisotropic stress is zero so that the Bardeen potentials are related by $ = — $ 
and V? is related to S, the comoving density contrast by [HI, l25j | 
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* = -3UrJ**' (2 - 2) 

where a is the scale factor, k the comoving wave number, H the Hubble parameter defined as H — a/a where the 
dot denotes differentiation with respect to coordinate time t, and we denote Fourier components by a subscript u k" . 
If inflation is driven by a light slowly rolling field <f> then, to leading order in the slow roll approximation, the field 
equation for field fluctuations, S<p, on flat slices lives in unperturbed spacetime. Defining ip = a5<f>, the Fourier 
components of the field perturbations obey 

d 2 ^(k,T) ( 2 



dr 2 
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fc 2 --U(k,r) = 0, (2.3) 



where r = —(aH) 1 is conformal time (it is assumed that inflation is exponential). The solution of this equation, 
with initial condition corresponding to the flat space-time mode function i/i(k, r) = exp (— ikr) / \/2k, is 

ip(k,T) = -4=(1 + ikr)-^— exp (-ikr) . (2.4) 
y/2k ikr 

The Fourier components of the comoving curvature perturbation 1Z are related to the field fluctuations on flat slices 
5(j) by 

n k =-^84> k . (2.5) 

For a definition of 1Z in terms of the metric perturbations and its relation to curvature perturbations defined on 
different hypersurfaces see for example Ref. [26|]. On sub- horizon scales therefore 

«, ^ 1 H 2 f k \ f ik\ . , 

The Bardeen potential * is related to 1Z by 

^9 k +<L*M* k = - il + w)1ikt (2.7) 



where w — p/p, with p and p being the pressure and the energy density, respectively. Under slow-roll conditions 
during inflation w ~ — 1, and therefore W is practically zero. 



We need to set the initial conditions for the radiation dominated epoch after the end of inflation. We assume for 
simplicity that the slow-roll conditions hold at the end inflation and that reheating is rapid. This is not necessarily the 
case; the decay and thermalization of the inflaton may be slow and in single field models slow-roll has to be violated at 
the end of inflation. Our assumptions are valid in, for instance, hybrid inflation models where the secondary waterfall 
field has a large mass so that the slow-roll conditions hold at the end of inflation, and the inflaton field rolls to the 
minimum of its effective potential at the end of the inflationary stage almost instantaneously [271 ] . 

Under these conditions, we can smoothly match the solutions of the Bardeen potential between inflation and 
radiation domination by requiring the intrinsic metric and the extrinsic curvature to be continuous on comoving 
hypersurfaces [28|, H^, [3(| • Then '5 during the radiation dominated epoch is given by [2l[ 

(r) = 3 ( 1 + ^fe( T c) [( x _ a; e ) cos (a; - a*) - (1 + xx c ) sin (x - x e )} , (2.8) 
2x^ 

where 

x = c s kr, (2.9) 

the subscript "e" denotes the end of inflation, and x e = c s kT e = c s k/a c H c . On super-horizon scales, k <C a c H e , and 
Eq. (|2.8j) reduces to the standard relation. On sub-horizon scales the Bardeen potential undergoes damped oscillations, 
reaching a maximum value ^(x*) at x — x+ during the first oscillation. The mass variance of is approximately 
given by a%(M) w Vy(M) so that 

ai(M, t) = AVn( f c ^ U x _ x \ C o S ( x - x c ) - (1 + xx c ) sin (x - x c )} 2 , (2.10) 
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where, for x c > c s , 



Vn{x c ) = An 



1 X e 



(2.11) 



defining A-r. = [(H% / '4> c ) /27r] 2 , the amplitude of the curvature perturbation power spectrum at the end of inflation on 
super-horizon scales. 

In Ref. [2l[, it was argued that Eq. (|2.8p is no longer valid when the amplitude of the curvature perturbation TZ 
becomes greater than one. Since TZ increases with (k/aH), there would then be a scale below which it appears that 
linear theory breaks down. At the end of inflation this scale is given by (k/a c H c ) ~ 1/TZnc, where TZnc is the 
amplitude of TZ at horizon crossing. However this is a conservative limit because well inside the horizon gravity should 
be less important, and therefore it is the scalar field perturbation becoming larger than unity which signals the break 
down of linear theory. An estimate of the critical wavenumber can be found by equating the energy density of the 
field perturbation and the background energy density. This gives (&/a c ) cr it ~ ^/M^H Cl which can be written in terms 
of the reheating temperature as 

The critical wavenumber decreases with increasing reheating temperature. The minimum value, which occurs for a 
reheating temperature of 10 16 GeV, is the same order of magnitude as the critical wavenumber corresponding to TZ ~ 1 
if 7?.hc ~ 0.01 — 0.1. For smaller reheat temperatures this criterion gives a larger critical wavenumber than TZ ~ 1. 
We discuss the consequences for the constraint on the amplitude of the curvature perturbation of these two limits in 

secnni 



B. PBH abundance 



In the standard super-horizon calculation of the abundance of PBHs one applies the Press-Schechter formalism, 
which is often used in large-scale structure studies (3lj . The quantity of interest (usually the density field, but in 
this application the Bardeen potential) is smoothed on a mass scale M, and the regions where the field perturbations 
exceed a certain threshold value are assumed to form PBHs with mass greater than M. The fraction of the Universe 



in PBHs with a mass similar to or larger than the smoothing scale 2 , is given by Press-Schechter theory, initially, as 
Op BH ,i(> M) ee F(M) = 2 / P(¥(M))dtt(M) = erfc ( J^Lj > ( 2 - 13 ) 



where cr\,(M) is the maximum mass variance of and the final equality follows from the fact that the smoothed 
Bardeen potential iff(M) has a Gaussian probability distribution. We have followed the usual Press-Schechter practise 
of including a factor of "2" in Eq. (|2.13|) . The collapse criterion in terms of this variable is \^\ > \^ c \ = 1/2, which is 
equivalent to a density contrast threshold at horizon crossing 5 C — 1/3. The mass fraction of PBHs per logarithmic 
mass interval is 

r dF(M) _ [2M\V C \ d<4 ( 9\ 



In the short wave-length limit, x c 3> 1, behaves like a damped sinusoidal function, and therefore the maximum 
occurs at x± — x c + (tt/2). In general there is no simple analytic solution for x+. The difference B(x c ) — (x* — x c ) 
converges very rapidly to the value n/2 inside the sound horizon. We used this fact in Ref. [2l[ to give an estimate 
of the mass fraction of PBHs. Here in order to extend the calculation to scales of the order of the horizon scale at 
the end of inflation, we calculate B(x c ) numerically. This calculation breaks down when the first maximum occurs 
outside the horizon, X* < c s . The PBH mass is related to x c by 

3* \kJ 3 h \kJ x c + B(x c ) C x 2 c (x c + B(x c ))' V ' 

where M c is the horizon mass at the end of inflation, 

with Trh the reheat temperature at the end of inflation and g c g the number of relativistic degrees of freedom; « 3 
and at high temperatures g c g ~ 100. 



III. CONSTRAINTS 



PBHs with Mpbh 5 x 10 14 g will not have evaporated by the present day (their lifetime is longer than the age 
of the Universe) and their density can not be too large. Traditionally 'not too large' was interpreted as S!pbh,o = 
PPBH,o/Ptot.o < lj where the subscript "0" denotes the present epoch. This constraint can be improved to Opbh o < 
^cdm,o ~ 0.47, where the numerical value is the 3-er upper limit on the CDM density extracted from WMAP [351 ] . 
however this factor of two improvement has an essentially negligible effect on the constraints on the primordial 
perturbations due to the exponential dependence of the PBH abundance on the size of the perturbations. 

Since ppbh oc a -3 and ptot oc a -3 (ptot °c a ) during matter (radiation) domination, f2pBH,o = ^PBH.eq and 

PPBH, q = r M^dM = r f^VdM, (3.1) 

J5xl0 14 g aM J5xl0 14 g aM V a cq / 

where the initial PBH number density dn;/dM is given by 

dn; Pi (M) dF(M) 



(3.2) 



Pi = ^ ff Jl 4 = P ra d, eq ^ =™ ^ h=H ■ (3.3) 



AM M dM 
The initial total (radiation) density is given by 



f „i T 4 / Scff W Ji \ _ Ptot.cq 9 C S \ J i 



On^en^l /-i^u.uq I cq I I rp I o I ct l / I 



2 We assume throughout that the PBH mass is equal to the smoothing scale mass, which in the standard, super-horizon case corresponds 
to the horizon mass Mjj • This is not strictly true, and numerical simulations indicate that the PBH mass depends on the size and shape 
of the perturbations 32, 33, 34]; however, this uncertainty is not important when applying PBH abundance constraints, due to their 
relatively weak mass dependence. 
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FIG. 1: The differential mass fraction of sub-horizon PBHs, eq. (|3.5[l . as a function of x e for A 1/2 = 0.02 and T RH = 1 GeV. 
The distribution is sharply peaked at small x D , demonstrating that the majority of PBHs form at essentially the same time. 

Combining Eqs. (13.2p . (j3.ip and (|3.3p . and using a oc g e g^ 3 T _1 , the present day mass fraction of sub-horizon PBHs 
becomes 



5xl0 14 g V a i 



a ea \ dF , , , 



(3.4) 



The upper mass limit of the integral in Eq. (|3.4p should correspond to a mass approximately given by the horizon 
mass at the end of inflation. However, the mass function decreases so rapidly for large masses that in practise the 
upper limit can be safely taken to infinity. Using Eqs. (|2.15|) and (|2.16|) the present day PBH mass fraction, can be 
expressed in terms of 



fipBHO = 7.2 x 10* 
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B(x e ) dx ( 
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where xq is the value of x e corresponding to M = 5 x 10 14 g which can be found from Eq. (|2.15p . For small reheating 
temperatures, Trh "C 10 5 GeV, we have B (xq) <C xq, and xq is given by 



r 3 

Xq 



46.2 x 



10 7 GeV 



(3.6) 



For larger reheat temperatures B(xq) can not be neglected, and therefore we calculate xq and hence B(xq) iteratively 
using Eq. (|2TT5|l . 

In Fig. [1] we plot the present day differential mass fraction of sub-horizon PBHs, df2pBH,o/da;c from eq. (|3.5p . as 

1 /2 

a function of x c for = 0.02 and Trh = 1 GeV. The PBH differential mass fraction is relatively sharply peaked 
with the majority of PBHs forming at essentially the same time. The present day PBH mass fraction, Eq. (|3.5p is 
therefore given, to a good approximation, by 



^pbh,o - 7.2 x 10 8 
for Xq 3> 1, and by 
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for xq < 1 where, in both cases, crj(xo) is the maximum value of the mass variance at x c — xq. For a range 
of reheating temperature values, Trh, we calculate the constraint on A-r, by first calculating xq, the value of x e 



corresponding to PBHs with M = 5 x 10 14 g (the lightest PBHs which will not have evaporated by the present day). 
We then find the mass variance of the Bardeen potential for this value of xq, cr^(xo), using Eq. ()2.8|) . Finally, we 
use either Eq. (|3.7p or (|3.8p . as appropriate, to calculate the constraint on An from the requirement fipsH.o < 0.47. 
The results are plotted in Fig. O As the reheat-temperature increases the constraint initially becomes tighter as the 
duration of the radiation dominated era, during which the fraction of the energy density in the form of PBHs grows 
as ppBH/ptot a, increases. As the reheat temperature increases the horizon mass at the end of inflation decreases 
and for 10 6 GeV < T RH < 10 8 GeV only a fraction of the sub-horizon PBHs have M > 5 x 10 14 g. For T RH > 10 s GeV, 
M c ~ 5 x 10 14 g and none of the sub-horizon PBHs are massive enough to last to the present day. Up to this point 
we have ignored the fact, as pointed out at the end of Sec. Ill Al that the calculation assumes that the linear theory 
holds. If we ignore PBHs formed on scales where 1Z > 1 at the end of inflation, we see in Fig. [5] that this weakens 
the constraints on A^[ 2 for Trh < 10 5 GeV, but the change is small (at most 20%). This is because the PBH mass 
function decreases rapidly for large x c . Taking the linear limit for the field perturbations in Eq. (|2.12[) the change in 
the constraint is negligible for the reheating temperatures for which our present day abundance constraint is applicable 
(T RH < 10 8 GeV). 

We now compare the sub-horizon constraints with the conventional super-horizon constraints. A full calculation 
would require the complete curvature perturbation power spectrum, or equivalently the specification of a concrete 
inflation model. In order to be as general as possible we simply consider the amplitude of the super-horizon curvature 
perturbation power spectrum, An, as a free parameter. On cosmological scales An = 2.4 x 1(U 9 [35(, therefore the 
production of an interesting (i.e. non-negligible) density of PBHs requires the curvature perturbations to be signifi- 
cantly larger on small scales. Assuming that the power spectrum grows monotonically with increasing wavenumbcr 
the abundance of super-horizon PBHs will be dominated by PBHs which form immediately after inflation at Trh- 
The present day abundance of super-horizon PBHs is then given by 

fWo(M > 5 x 10 14 g) = 7.2. x 10 8 -J|erfc ( Sc ) . (3.9) 

bev V v2cT(5(Jrh)/ 

The Fourier components of the density contrast and the curvature perturbation are related by [36j ] 



2(1 +w) ( k 
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^=t^UJ n » (3 - 10) 

so that, at horizon crossing, V s = (4/9) 2 P K = (4/9) 2 „4 TC . Using a s (M) w V S (M) again, then a s {M) = (4/9) 2 An- 
The threshold density has historically been taken as S c — w, where w = p/p is the equation of state [2|, so that 
for formation during radiation domination (which is usually the case of interest) 5 C = 1/3. More recently Green et 
al. [37l ] carried out a new calculation using peaks theory and a formation criterion derived by Shibata and Sasaki [32j 
which refers to the central value of the metric perturbation. They found that the standard Press-Schechter calculation 
agrees with the peaks theory calculation if the density threshold contrast is in the range 0.3 to 0.5, while numerical 
simulations by Musco et al. found a threshold S c = 0.45 [Hj]. We therefore consider two density thresholds, S c = 1/3 
and 0.45. 

The constraints on the curvature perturbation at the end of inflation, as a function of reheat temperature, are 
plotted in Fig. [H The super-horizon constraints are a factor of ^3 — 4 weaker than the sub-horizon constraints, 
and the factor of ~ 1.4 difference in the density contrast thresholds we consider gives, as would naively be expected, 
a similar difference in the constraints (with the larger threshold leading to a weaker constraint). The super- horizon 
constraint stops at Trh ~ 10 8 GeV as then M c < 5 x 10 14 g and the PBHs which form immediately after inflation have 
evaporated by the present day. PBHs with M > 5 x 10 14 g could form later, at T < Trh, though their abundance 
would depend on the precise scale dependence of the power spectrum. 



IV. DISCUSSION 



CMB and LSS observations probe primordial perturbations on scales of the order 1 Mpc to 3000 Mpc that typically 
leave the horizon roughly 50-55 e-foldings before the end of inflation. The amplitude and scale dependence of the per- 
turbations on these scales are now determined fairly accurately through observations, however these scales correspond 
to only a small region of the inflationary potential. On the other hand PBHs constrain the primordial fluctuations 
over a large range of far smaller scales, and hence provide additional and complementary constraints on the dynamics 
of inflation over a wider range of energy scales. These constraints are relevant for models where the perturbations 
have significant scale-dependence, for instance models where there are multiple stages of inflation or inflation ends 
via a phase transition. 
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FIG. 2: The constraint on the amplitude of the curvature perturbation power spectrum at the end of inflation from the present 
day abundance of PBHs as a function of reheat temperature Trh- The solid and dotted lines are for sub-horizon PBHs. The 
solid line assumes that the linear theory calculation is still valid on scales on which the comoving curvature perturbation 
becomes larger than one before the end of inflation, while the dotted line ignores PBH formation on these scales (and hence 
provides a conservative evaluation of the constraint). The long and short dashed lines are for super-horizon PBHs using a 
density contrast threshold S c = 1/3 and 0.45 respectively. 

Observational constraints on the abundance of PBHs are usually carried out for scales which exit the horizon 
during inflation (with PBH formation occurring after horizon re-entry during the subsequent radiation dominated 
era). Lyth et al. 2l| recently studied the formation of PBHs on smaller scales on which the perturbations do not exit 
the horizon during inflation, and therefore never become classical. We have extended this work by calculating the 
constraints on the amplitude of the primordial curvature perturbation power spectrum at the end of inflation, An, 
from the present day abundance of these sub-horizon PBHs, as a function of the reheat temperature after inflation. 
For Trh < 10 6 GeV most of the sub-horizon PBHs have M > 5 X 10 14 g and hence have lifetimes longer than the 
age of the Universe. The requirement that their present-day density does not exceed the upper limit on the cold 
dark matter density leads to the constraint An < 0.25 — 0.3 with a weak dependence on the reheat-temperature (the 
higher the reheat-temperature, the longer the duration of the radiation dominated epoch during which the fraction 
of the energy density in PBHs grows). As the reheat temperature increases the horizon mass at the end of inflation 
decreases. For 10 6 GeV < Trh < 10 8 GeV the fraction of the sub-horizon PBHs which have M > 5 x 10 14 g decreases 
rapidly with increasing Trh and hence the constraint on An is weakened. For Trh > 10 8 GeV, M ~ 5 x 10 14 g so 
that none of the sub-horizon PBHs are massive enough to last to the present day and hence there is no constraint on 
An from the present day abundance of sub-horizon PBHs. 

We also calculated the constraint from the present day abundance of PBHs which form on larger scales, that do exit 
the horizon during inflation. This is a factor of ~ 3 — 4 weaker than the constraint from sub-horizon scales. For higher 
reheat temperatures there will be constraints from the effects of the PBH evaporation products on the successful 
predictions of nucleosynthesis @, H i, 0, H, S El EL H? ]■ 

These constraints are however model dependent, 
depending on the PBH mass function and hence the scale-dependence of the primordial power spectrum, which is 
beyond the scope of this paper. 

There may also be constraints from the present day density of relics which may be the end point of PBH evaporation 
and the emission of long-lived massive particles. The applicability, or otherwise, of these constraints depends on the 
details of physics beyond the standard model of particle physics and are hence also beyond the remit of this paper. 
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